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Abstract — This paper explores using RSS measure- 
ments on many links in a wireless network to esti- 
mate the breathing rate of a person, and the location 
where the breathing is occurring, in a home, while the 
person is sitting, laying down, standing, or sleeping. 
The main challenge in breathing rate estimation is 
that motion interference, i.e., movements other than 
a person's breathing, generally cause larger changes 
in RSS than inhalation and exhalation. We develop 
a method to estimate breathing rate despite motion 
interference, and demonstrate its performance during 
multiple short (3-7 minute) tests and during a longer 
66 minute test. Further, for the same experiments, 
we show the location of the breathing person can be 
estimated, to within about 2 m average error in a 56 
m^ apartment. Being able to locate a breathing person 
who is not otherwise moving, without calibration, is 
important for applications in search and rescue, health 
care, and security. 



I. Introduction 

This paper explores using standard wireless devices 
which measure only received signal strength (RSS) to 
monitor and localize the breathing of a person, who 
does not wear any device, in a residential environment. 
In comparison, existing breathing monitoring methods 
either require sensors in contact with the body, or require 
specialized and expensive radar devices [1], [2]. The core 
of the idea of non-contact breathing monitoring via RSS 
is that some links' RSS values, on a particular channel, 
are highly sensitive to motion near the link line (the 
line between transmitter and receiver), and that those 
links can be used, in the absence of other motion, to 
estimate the person's breathing rate, and in addition, to 
estimate that a breathing person is near that link line. 
From many links' RSS data, and from measurements on 
multiple frequency channels, we can estimate breathing 
rate and perform breathing localization. One of the core 
challenges is that these same links also experience rapid 
change due to other movement, for example, if a person's 
limb moves. Our methods are robust to these normal 
occasional movements because it identifies when a sudden 
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movement occurs, subtracts its effects out, and can still 
then estimate breathing rate in periods containing that 
motion. Our work is the first, to our knowledge, to use RSS 
measurements to locate where that breathing is occurring. 

A significant body of related research has shown that 
people moving people can be located in a building using 
a static wireless network that measures links' RSS values 
[3], [4], [5], [6], [7], [8], [9]. This can be done through walls 
[To] , [H] , [12] , using a variety of statistics of the measured 
RSS. However, state-of-the-art RSS-based non-cooperative 
localization methods require either: 1) motion, that is, a 
person to be moving at least once during the course of 
the time in which measurements are collected [10], [13j . 
14 , or 2) empty-area calibration, measurements of the 
RSS from the period of time when the area was cleared of 
any people [5], [5]. Even systems sophisticated enough to 
adaptively learn to distinguish the statistics of RSS during 
crossing vs. no crossing require some periods of both to 
occur jl2j . In an search and rescue operation or a building 
fire, people may be unconscious and thus motionless by 
the time emergency responders arrive, and thus a non- 
cooperative user localization system may not be able to 
locate them. Although some work on robotic RF sensor 
networks has demonstrated the ability to image concrete 
structures |15j . work has not yet shown that a motionless 
person could be located without prior empty calibration. 
In this paper, we show a method that would locate a 
stationary but breathing person without an empty-area 
calibration, although we note that we do not test it 
through walls or in rubble. 

As another application, using standard wireless net- 
works for contact-free breathing monitoring may be useful 
for health care and wellbeing. One may be able to use 
a wireless network deployed in a bedroom as a baby 
breathing monitor, or as part of an in-home sleep apnea 
diagnostic system. Elder care systems which track and 
monitor an elderly person's location and activities for 
abnormal changes could use breathing rates as additional 
context information. As of yet, we have not developed algo- 
rithms to identify or distinguish the breathing rates of two 
people in the same room. However, in these applications, 
breathing monitoring is most critical when the person is 
alone, and thus no other person is present to notice their 
apnea or their stoppage of breathing. 

Finally, in general, breathing rate provides a measure 
of a person's context, and context-aware computing sys- 
tems could benefit. Perhaps "smart" homes could also be 
"empathetic" and respond appropriately when your at-rest 
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breathing rate is abnormally high. A scary movie could 
adapt to the viewer's measured level of stress. 

In the area of contact-free breathing monitoring, most 
systems use specialized radars. Doppler radar devices can 
be used p], but slow and small breathing movements 
mean that the Doppler shift is very low unless the center 
frequency is very high. As a result, commercial Doppler 
radar breathing monitoring products are expensive and 
out of reach for many applications. For example, the Kai 
Medical "Continuous" system [TB], although not yet FDA 
approved, is said to be priced at $2000. Ultra wideband 
impulse radars can estimate breathing as well as heart rate 
[U . Radar propagation losses are of order in free space 
[17] , as opposed to cP for transmission TS] , and thus radar 
devices either have limited range or high transmit power. 
Regarding the latter, UWB devices are severely limited in 
transmit power by regulation, and thus compliant devices 
have a low range. We believe that use of a standard 
wireless network offers a low-cost alternative for residential 
breathing monitoring. 

Our initial work was presented in |19| , which introduced 
RSS-based breathing monitoring, presenting algorithms 
for estimating breathing rate and detecting the absence of 
breathing. The methods were experimentally verified on a 
subject in a hospital bed, using transceivers placed within 
centimeters of the subject's body. The report shows that 
links in a deep fade are the ones most likely to measure 
changes when the person breathes. Our breathing rate 
estimator is compared to the method of [T9' in Section 

imi 

This paper makes several advancements: 

1) This paper is the first to demonstrate that breath- 
ing can be monitored using wireless networking 
devices meters, as opposed to centimeters, away 
from the person. In addition, we use measurements 
on multiple frequency channels, and quantify the 
improvement possible compared to single-channel 
RSS measurements in Section IV-Bl Our testbed and 
experimental setups are described in Section HIl 

2) This paper presents an algorithm to estimate and 
avoid motion interference, and as a result, estimate 
breathing rate even during windows of time in which 
there is motion in the environment. The subject 
in |19) was completely stationary, except for an 
expanding and contracting chest. When the person 
moves, for example, coughs, scratches an itch, moves 
an arm or leg, or changes position in bed, the method 
of [19] becomes unable to estimate breathing rate 
for a period of time {e.g., 30 s). Fundamentally, 
this is because the links most sensitive to chest 
motion are also those most sensitive to other motion. 
Our method shows dramatic breathing rate estima- 
tion performance improvements. Using 30 seconds 
of data, our method's breathing rate estimates have 
average error of 1.0 bpm, compared to 1.7 bpm for 
[T5] . Our breathing rate estimation algorithm and 
experimental results are presented in Section IIIIl 

3) Third, the paper is the first to develop and test 



an algorithm to locate a breathing but otherwise 
stationary person. Without any empty-room calibra- 
tion or training period, we estimate an image map 
of breathing intensity in the deployment area, and 
locate the breathing from the highest point in the 
image. We see average localization errors of about 
2 m. Our localization algorithm and experimental 
results are presented in Section ITVl 

Finally, we conclude and discuss future research directions 
in Section IVll 

II. Experiments 

In this section, we describe the testbed network, the 
residential environment in which the tests are conducted, 
and the testing procedure. 

A . Network 

In general, we suppose there is a network of S 
transceiver nodes which can operate on one of C frequency 
channels. Link I has transmitter (TX) Ti, receiver (RX) Ri, 
and frequency channel Fi, i.e., each TX / RX / channel 
combination is a different logical link. In our experiments, 
L = CS{S — 1) because all links are connected, but full 
connectivity is not necessary for breathing monitoring. 

We assume that each link's RSS is measured at regular 
intervals with period T by sending and receiving a packet. 
Note that newborns have the highest breathing rates, near 
37 breaths per minute (bpm), or 0.62 Hz, and adults at 
rest breathe near 10-14 bpm, or 0.23 Hz [20]. To sample 
at the Nyquist rate, sufficient for a newborn monitoring 
application, our sampling rate would need to be more than 
1.24 Hz. We denote the nth RSS on link / as ri[n], in dB 
units. Future work could extend the presented algorithms 
for irregularly sampled data. 

In our experiments, we use the TI CC2531 dongle node 
21|, an IEEE 802.15.4 compliant radio. The CC2531 
transmits with transmit power 4.5 dBm in the 2.4 GHz 
ISM band, and can select from one of 16 frequency chan- 
nels. 

A multi-channel TDMA protocol is used in which each 
node transmits in sequence, and after all nodes' slots, 
nodes switch synchronously to the next frequency channel 
[35]. A node's packet transmission includes its node ID 
and the RSS of the most recent packet received from each 
other node. A sink node that overhears all packets stores 
the data for post-processing. 

B. Apartment Experiment 

In our first experiment, the apartment experiment, we 
deploy 5* = 33 nodes, each at a height of 0.9 meters. Our 
protocol uses C = 4, specifically, 802.15.4 channel numbers 
15, 20, 25, and 26, selected to overlap the least with WiFi. 
Each link is measured each T = 428 ms. This is a sampling 
rate of 2.3 Hz. 

The testbed network is deployed in the apartment shown 
in Figure |T] Node locations are indicated with circles on 
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Fig. 1. Floor plan of 7 X 8 m apartment, deployed nodes (•), and 
dimensions in cm. 



the map. The apartment is on the third floor of a fully- 
occupied six story building in downtown Salt Lake City, 
Utah. It is an actual person's residence, and is furnished, as 
can be seen from the photo of the living and kitchen areas 
in Figure [5] Tests are conducted on a Saturday between 
noon and 1:00 pm. During the tests, other people are 
sometimes heard in the hallway outside of the apartment 
and in neighboring apartments. 




Fig. 2. Living and kitchen areas in apartment, with nodes seen taped 
to the wall at height 0.9 m. 

As a large apartment building in a metropolitan area, 
there are many sources of 2.4 GHz interference. The 
apartment has its own wireless router, and more than ten 
strong WiFi access points are observed, and although the 
channels are chosen to have the least overlap with WiFi, 
they are not perfectly orthogonal. Microwave ovens and 
cordless phones also exist in the nearby environment. 
Tests: Five tests are conducted. In each, the experimenter 
stays in a single location for 3-7 minutes and maintains 



maintains a constant breathing rate of 10 bpm using 
a metronome. The experimenter is instructed to sit or 
stand still at the location, but to move as required to be 
comfortable, i.e., move an arm or leg or change position. 
The experimenter made some movement, on average, twice 
per minute. In the tests, the experimenter is located: (1) 
standing in the kitchen; (2) sitting by the dining table; (3) 
sitting on the sofa; (4) laying in bed; and (5) sitting on the 
toilet. 

C. Nap Experiment 
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Fig. 3. Master bedroom floor plan, for "nap" experiment. Where two 
nodes are at the same (x, y) coordinate on this 2D map, the higher 
node number is listed first. 

Our second experiment, the nap experiment, is con- 
ducted in the master bedroom on the second floor of a 
house, as shown in Figure [31 S" = 12 nodes are deployed in 
a 5.3 m by 5.3 m area, at power outlet locations or in power 
strips in the room. We placed two nodes per power outlet, 
one on the floor and the other at a height of about 0.48- 
1.65 ni. The height of the high nodes varied because they 
are placed on top of furniture (bedside table, dresser, table, 
etc.). In this experiment, C = 5 channels are measured, 
specifically, IEEE 802.15.4 channels 11, 15, 18, 22 and 26. 
Here, the sampling period was 179.6 ms, for a sampling 
rate of 5.56 Hz. 

There is also WiFi interference at this location, with 
one wireless router on the first floor below the bedroom. 
A laptop at location (2m, 0.5m) records data and sends 
it over WiFi to a server, ensuring some level of WiFi 
interference to the breathing monitoring system. 

A single 66-minute long test is conducted. The experi- 
menter starts collecting data on the laptop, walks to and 
climbs into the bed, and takes a nap. In contrast to the 
previous experiments, no ground truth is known for the 
experimenter's breathing rate, as no separate breathing 
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monitoring sensor is attached to the person, and the 
person cannot maintain a constant breathing rate while 
sleeping. After time 65 minutes, the person wakes up, gets 
out of bed and stops the data recording. The person moved 
at least once during sleep, since it was noted that he did 
not wake up in the same position that he went to sleep. 

III. Breathing Rate Estimation 

A. Basic Method 

In [12], an RSS-based breathing rate estimator is pro- 
posed which, in short, calculates the power spectral density 
(PSD) over each link using the most recent N samples, 
sums the PSD over L links, and estimates the breathing 
rate as the frequency at the maximum of the sum PSD. 
Mathematically, 
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f = argmax 



E 

1=1 



yi [n\e 



n=i-N+l 



(1) 



where fmin and fmax are the minimum and maximum 
breathing rates which the system is set to allow, i is the 
current time, j — ^/—l, and yi[n] is the change in RSS 
compared to the mean, defined as 



yi[7i\ = ri[n\ - ri, 



(2) 



where n is the average RSS on link /. We call yi [n] the RSS 
signal because the mean RSS carries no information about 
breathing, and in fact, a strong DC component overwhelms 
any other frequency contained in the data. In |19| . DC 
removal is performed using a 7-tap IIR filter. However, we 
find that setting ri — X]n=i-Af+i i-e., the window 

average, has superior performance over a wide variety 
of experiments we have conducted. When comparing our 
method, to be fair to the spirit of the method in [19], we 
use this window average in ([2|). 

B. Breakpoint Method 

The challenge faced by the RSS-based breathing mon- 
itoring in general, and the basic method in particular, 
is that the links which best measure the person's chest 
movement are also those that best measure other motion in 
the environment. This other motion, which we term "mo- 
tion interference", can hide the breathing-induced changes. 
From our observation, motion interference and breathing 
cause notably different types of changes. While breathing 
causes slow, periodic RSS changes, motion interference 
typically causes sudden changes. For example, a person 
who rolls over in bed, or moves a foot, causes fast RSS 
changes during their movement. When the person stops 
moving, a link's RSS typically settles to an average value 
different than before the movement. If such a shift occurs 
during the iV-sample window, subtracting the mean as 
in ([2|) will not suffice to remove the strong low-frequency 
component in the RSS signal. 

Instead, we propose to detect each time index during 
which a sudden RSS change occurs. We call these detected 
time indices "breakpoints". At each breakpoint, the system 



Parameter 


Value 


Window Duration, A' 


70 (30 sec) 


Breathing Rate Min., fmin 


0.1 Hz (6 bpm) 


Breathing Rate Max., fmax 


0.4 Hz (24 bpm) 


Nodes, S 


33 


Channels, C 


4 


Sampling Period, T 


0.428 sec 


T-test Window, Q 


14 (5 sec) 


RMS T-score Threshold, 7 


0.8 



TABLE I 

Breathing Rate Estimator and Experiment Parameters. 



is allowed to forget the previous link RSS averages and 
calculate new averages. By doing so, we are able to perform 
more accurate mean removal, and thus make breathing 
rate estimation more robust to motion interference. 

The T-test is an appropriate statistical detection test to 
perform this task. The T-test is a generalized likelihood 
ratio test (GLRT) of the change in mean between two 
groups of data [23], in our case, one group of RSS signal 
samples prior to the timc-under-test (the sample index at 
which we want to detect a change in mean), and another 
same-duration group of data after the time-under-test. 
Specifically, the group t-score is. 



ri [n\ - ri [n\ 



max 



A)Iq}' 



(3) 



for each link I and time t, where n[n] and r;[n] are the 
average of the Q samples of ri[m] before and after time 
n, respectively, crf[n] and ijf[n] are the sample variances 
of r; [m] before and after time n, respectively, and e > 
is used to prevent division- by-zero. We use e = 0.5 in our 
experiments. 

Next, at each time n we compute the root-mean-squared 
(RMS) average of Ti[n] over all links Z, 
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TRMs[n] = I j^^T^W] 

\ 1=1 



(4) 



Time n is a breakpoint index if TRMsi^] — 7j where 7 is 
a predetermined threshold. For any time window, we also 
consider the starting time index and ending time index as 
breakpoints. At the current time i, with a window length 
N , the starting and ending time indices are i — iV -I- 1 and 
i. 

Breakpoints are used to remove the mean from the raw 
RSS as follows: 



ri[n\ - ri[n\, 



(5) 



where f; [n] is the average of ri [m] for all m such that bp < 
m < bf, where bp is latest breakpoint before or at time 
index n, and 6/ is the earliest breakpoint after time n. 
Essentially, the mean model f; [n] is piecewise constant, 
with transition times at each breakpoint. After performing 
mean removal in ([S]), we use our new RSS signal yi[n] in 
([1} to perform breathing rate estimation. 
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Fig. 4. The mean-removed RSS signal yi [n] for four links that 
best show breathing during 'sofa' experiment, using (top) basic and 
(bottom) breakpoint methods, with estimated breakpoints (:). 

C. Apartment Experiment Results 

To compare the basic and breakpoint methods, we 
describe an example taken from a 30-second window 
during the 'sofa' experiment. Methods are tested using 
parameters given in Table ID Figure |4] shows the mean- 
removed RSS signal yi[n\ for the basic and breakpoint 
methods, for four links I that are the four best in terms 
of signal amplitude at the true breathing frequency. In 
the basic method, yi [n] drifts upward, particularly around 
n = 213. Many links (including those not shown) exhibit 
a shift in mean here, as well as near n — 175. The 
breakpoint method, using a threshold 7 = 0.8, calculates 
a high RMS t-score near these two times, and detects 
several time indices near 213 and 175 as breakpoints. By 
removing the mean in between each pair of breakpoints, 
the signal yi[n] no longer has a drifting mean in the 
breakpoint method. Note that some ranges are reduced 
to zero because breakpoints are dense, and removing the 
average from a single sample results in a zero value. 

Figure [S] shows that the breakpoint method enables 
correct breathing rate estimation when the basic method 
is unable. The basic method's PSD is shown to be highest 
at the lowest breathing rate, 0.10 Hz, and thus / = 0.10 
Hz. The PSD for the breakpoint method is reduced at all 
frequencies, but has a clear maximum at 0.166 Hz, very 
close to the true 0.167 Hz breathing rate. 

Next we test all 30 second windows for all test locations, 
moving the window 5 seconds each time. We first report 
the fraction of estimates / which are acceptable, which 
we arbitrarily define as being within 3.0 bpm of correct. 
Second, we report the average frequency error over all 
windows. 

Results are shown in Table HIl Over all tests, the basic 
method obtains an acceptable rate estimate 61% of the 
time. The breakpoint method achieves 81% acceptable es- 
timates, approximately cutting in half the number of 'un- 



Fig. 5. Normalized average PSD from ((TJ using basic and breakpoint 
methods and data from Figure|4] estimated breathing rate (• and □), 
and true breathing rate (:). 





% Acceptable 


Avg. 


Err. (bpm) 


Test Loc 


Basic 


Breakpt 


Basic 


Breakpt 


Sofa 


71% 


94% 


1.23 


0.348 


Table 


81% 


89% 


1.14 


0.831 


Kitchen 


58% 


91% 


1.81 


0.859 


Bathroom 


71% 


97% 


1.29 


0.339 


Bed 


25% 


33% 


2.95 


2.63 


Average 


61% 


81% 


1.69 


1.00 



TABLE II 

Basic vs. breakpoint breathing rate estimation: percent of 
estimates "acceptable'', and average error. 



acceptable' rate estimates. With the breakpoint method, 
four of five tests have acceptable rates above 89%, whereas 
the basic method had none. The breakpoint method re- 
duces the average breathing rate error to 1.0 bpm, a 41% 
reduction from the basic method. Primarily, the reduction 
in average error comes from eliminating large errors - both 
methods typically have less than 0.3 bpm error when the 
rate estimate is "acceptable". For reference, an end-tidal 
CO2 meter, the gold standard breathing rate monitor in 
hospitals [23], is accurate to ±1 bpm. 

The bed test has particularly poor performance. We 
note that if we had used the breakpoint method with 
7 = 0.5 (instead of 0.8), the bed test would have had 
73% acceptable estimates (vs. 33%), however, the table 
test would fare worse with the lower threshold. Future 
work should address adaptive methods to decide upon 
breakpoints. 

D. Nap Experiment Results 

The data from the nap experiment is processed with the 
breakpoint method, using the same parameters as given in 
Table H 

As discussed in Section III-C[ no ground truth is avail- 
able, however, there are features of the breathing rate 
estimates which match our understanding of the human 
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Fig. 6. Breakpoint metliod nap experiment breathing rate estimates 
based on 30 s data window (■), 90 s median ( ), and breakpoints 

( ) which indicate periods when movement is detected. 



the number of rate estimates at or very near 6 bpm as 
a measure of the abihty of the method to remove normal 
(non-breathing) motion from the breathing-induced signaL 
In this experiment, using the breakpoint method, there 
are fourteen breathing rate estimates below 9 bpm, five of 
which are during identified periods of motion. Thus out 
of 793 window / estimates, only nine (1.1%) are very low 
and are during periods when no motion is detected. 

In comparison, we show results for the basic mean- 
removal method in Figure[71 When using the basic method, 
two metrics are worse. First, there are 48 estimates with 
rate below 9 bpm. Since there are several of these low 
estimates in a row, the median filter also falls below 9 
bpm during periods that we believe (from the breakpoint 
method results) contain movement. In contrast, the me- 
dian filter in the breakpoint method never falls below 
9 bpm. Not only are the basic method rate estimates 
presumably incorrect, in addition, there is no "motion" 
label provided to indicate to the application that there was 
likely motion in that period. As a result, there is no way 
to automatically discard the poor rate estimates returned 
by the basic method. 



body during sleep. It has been shown that breath rate 
decreases during the transition from the awake state to 
stage two of sleep [25]. This is primarily because there 
is a downward trend in energy expenditure during the 
onset of sleep [25], [26]. The breath rate in Experiment 
2 decreases approximately 20% in the first 20 minutes, 
which correlates with published data showing a measured 
metabolic rate decrease of approximately 20% during the 
first 20 minutes of the onset of sleep |26j. During sleep 
stages 2, 3/4 and REM (rapid eye movement), the venti- 
lation and metabolic rate are known to become relatively 
stable [IS] , [17] , which is consistent with what we observed 
in our nap experiment estimates. However, moving during 
sleep increases energy expenditure and thus could tem- 
porarily increase breathing rate [25], as we see clearly at 
time 38 minutes, just after a period in which motion is 
detected by the breakpoint method. At the end of the nap 
experiment, the breath rate increased as the person awoke 
without an alarm. This is consistent with published data 
showing increased movement and metabolic rate during 
the transition to wakefulness [261 . 

In a period containing motion that is not detected and 
removed by the breathing rate estimation algorithm, the 
rate estimate tends to rail to /mm- This is because a 
single step function, when the step is large compared to 
other signals, would be estimated to be evidence of one 
cycle in the measurement period. For a 30 second window, 
this would be a rate of 2 bpm. Since fmin = 6 bpm 
in these experiments, the estimator will return a value 
of 6 bpm in this case. Step functions are not removed 
when using the basic algorithm, and in some instances, the 
breakpoint method may not detect a step function if the 
amplitude of the change is small compared to the variances 
before and after the jump. Although we are not certain 
of the sleeping person's true breathing rate, we quantify 




20 30 40 50 

Time Elapsed (minutes) 

Fig. 7. Basic method nap experiment breathing rate estimates based 
on 30 s data window (•), and 90 s median ( ). 



IV. Breathing Localization 

In this section, we describe a method to compute a 
"map" of where the breathing, at the estimated rate, is 
occurring. From the map we estimate the coordinate of 
the breathing experimenter within the apartment. 

Variance-based radio tomographic imaging (VRTI) uses 
the variance of RSS on links in a wireless network to 
identify the location of a moving person in a building 
[To] . However, a breathing but stationary person does 
not show up in a VRTI image, as the image is more 
affected by noise than by the small changes caused by 
breathing. In this section, we adapt the approach from 
10 for breathing mapping and localization, and we show 
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Parameter 


Value 


Pixel Width, Sp 


0.2 in 


Pixel Variance, cr^ 


2 


Correlation Distance, S 


2 m 


Ellipse Size Parameter, Ae 


1 m 


TABLE III 



Breathing mapping and localization parameters 





Avg. Loc. Error (m) 


Test Loc 


Basic 


Breakpoint 


Sofa 


1.5 


1.7 


Table 


2.1 


2.1 


Kitchen 


2.9 


3.6 


Bathroom 


1.5 


1.9 


Bed 


2.6 


2.7 


Average 


2.1 


2.4 



TABLE IV 

Average location errors for basic and breakpoint methods 



using the experimental data that an approximate location 
of the breathing person can be determined. 



A. Method 

We use the PSD at the estimated breathing rate for each 
link as an input to our method. For link I, at time i, the 
PSD of link I at / is, 



vi 



J2 yiHe-^''''^- 

n=i-N+l 



(6) 



The vector of all link values \s v = [vi, . . . ,vl\'^ . 

We want to estimate an image of breathing amplitude 
vs. space, denoted vector x, where Xk represents the 
quantity of breathing energy coming from pixel k. As in 
[lOj . we assume that w is a linear combination of x via a 
weighting matrix W plus noise 77, 



V = Wx + rj 

where W is defined as having (Z, fc) element given by, 



I|ZT,-Pfc|| + ||ZR,-Pfc|| 



Pi ' 

0, o.w. 



< 1 



(7) 



(8) 



where zt, and zn^ are the coordinates of TX and RX for 
link I, respectively, pk is the coordinate of pixel fc, Ae is the 
ellipse size parameter, and P; is the number of non-zero 
weights for link I. Essentially, Pi normalizes the weight so 
that the total weight of each link is 1. 

Solving for a; in ([71) given measurement v is ill-posed. 
We assume the image vector has covariance matrix G, 
with (fc,m) element Gk,m = cr^e""^*"^""/'', where is 
the variance of any element of x, and d is the correlation 
distance. The regularized least squares solution for x is 
thus. 



IIu, where 11 



(9) 



Note n must be computed only once. The real-time com- 
putation of the image requires only one matrix multiply, 
of O {LP) multiplies and adds, where P is the number of 
pixels. 

From the image, the coordinate of the pixel with max- 
imum value in x is used as the location estimate for 
the (single) breathing person. Future work must address 
localization in the multi-person case. 



B. Apartment Experiment Results 

We implement the breathing mapping and localization 
of (0) using the parameters given in Table Hill and evaluate 
its performance using the apartment experiment. Using 
the same example window as that shown in Figure |4l we 
first draw link lines corresponding to the few links I with 
highest VI in Figure IHlJb). These include the four plotted 
in Figure SI The image x for the same data is shown in 
Figure |Sl[c). For this example, the location estimate has 
0.82 m error. 

We summarize the performance of the localization es- 
timator across all windows and tests in the apartment 
experiment in Table HVl Using the basic method for mean 
removal, the average localization error over all tests is 
2.1 m. In a 7 by 8 meter area, this is a somewhat 
coarse estimate, however, it is typically sufficient to tell 
which room the person is in. Counterintuitively, for the 
breakpoint method, the localization estimates degrade 
somewhat to 2.4 m. From examination of the RSS data 
and images produced using the basic method for mean 
removal, we can see that when windows include motion 
interference caused by the experimenter, the breathing 
map X can be as good or better than without the motion 
interference. This is because the experimenter is causing 
the motion, which then leaks into the RSS signal on links 
near the experimenter, and helps to increase the image 
value near the experimenter's location. The breakpoint 
method removes much of this motion interference, which 
then prevents it from increasing the image value near the 
correct location. 

We notice some times when the image suddenly shifts 
and shows a bright spot at the top of the apartment, 
where there is a hallway outside of the apartment wall. 
The experimenter had heard people walking through the 
hallway during some tests. However, we did not record 
the ground-truth times of these events and thus future 
controlled experiments are suggested. 

Finally, we note computation time is about 0.2 s per 
30 s window, for both rate estimation and localization, in 
Python on a two-core 2.0 GHz processor - thus real time 
monitoring is very possible. 

C. Nap Experiment Results 

Finally, we test the nap experiment data for localization 
performance. We use the same algorithm and parameters. 
There is only one actual coordinate, which we measure 
to be the position of the chest when the person started 
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Fig. 8. (a) lines for the top ten links by Vi, and (d) breathing image x (white=high), with true (x) and estimated (+) locations. 



napping. There are 792 estimates, which are shown in 
Figure [9] Since we use the center of the maximum pixel of 
X as the location estimate, many of the location estimates 
exactly overlap. We also show the average of all location 
estimates in Figure [9l This average is about 1.0 m away 
from the true chest position. Considering the estimates 
individually, the RMS localization error is 1.16 m. In 
general, performance is better than the average of the 
apartment experiment, but the size of the area is also 
smaller as well, about 1/3 of the area of the apartment. 



V. Performance vs. Size 

In this section, we explore the performance of breathing 
monitoring when a smaller number of sensors, or a smaller 
number of channels, or both, are used in the network. 
From a system design perspective, it would be beneficial to 
deploy a smaller number of nodes, and for those nodes to 
make as few measurements as possible. If fewer channels 
were required to be measured, then one might transmit 
and receive less, thus enabling lower energy consumption, 
and less interference to other devices. 
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Fig. 9. Location estimates for each 30 second window (of which there 
are 792) during the nap experiment, along with an average over all 
windows, and the actual coordinate of the center of the person's chest 
during the nap. 



A. Number of Nodes 

In this section, we show the performance of breathing 
monitoring when we only use a subset of nodes, to test 
using experimental data how the performance changes as 
a function of S. For this, we use the data from the nap 
experiment. We would like to know, how many nodes are 
required for acceptable performance? 

Since we do not have a "ground truth" breathing rate 
for the nap experiment, we evaluate the performance of 
breathing rate estimation as follows. We assume that 
short-term variation in breathing rate estimates are char- 
acteristic of a less-accurate estimator. As such, we report 
the RMS difference between a current breathing rate and 
the short-term median, where the median is taken from all 
breathing rates estimated in a 90-second window centered 
at the current time, which we refer to as the "RMS- 
median". 

First, we test using random subsets of the twelve nodes. 
For each number of included nodes S < 12, we run 30 
random trials. We show the RMS-median value for each 
trial, as well as the average vs. S, in Figure [101 The RMS- 
median value, on average, decreases with S; however, the 
lowest possible RMS-median value seems to have a shallow 
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Subset 


RMS-median 


of Nodes 


(bpm) 


{0,2,4,6} 


0.72 


{1,3,5,7} 


1.96 


{0,2,4,6,8} 


0.84 


{0, 1,2,3} 


3.57 


{0,1,. ..,5} 


1.31 


{0,1,..., 7} 


0.68 


{0,1,..., 9} 


0.73 


{0,1,..., 11} 


0.91 



TABLE V 

RMS-MEDIAN USING NODE SUBSETS FOR NODES NUMBERED AS IN 

Figure m 



minimum at a value of S — 6, and every tested subset size 
had some trial with RMS-median value below 1.0. It would 
seem that low number of nodes is fine if you select the right 
nodes. 
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due to the measurement of multiple channels. Here, we 
test the performance during the nap experiment when only 
links from a subset of the channels are used. We test each 
possible (non-empty) subset of the five measured channels, 
and plot the RMS-median performance in Figure [TTJ First, 
using all 5* = 12 nodes, we see average performance 
improves as the number of utilized channels increases, 
although slowly. At C = 1, the average RMS-median 
is 1.56 bpm, while when using all C = 5 channels, the 
value is 0.91, a 41% reduction. The result is much more 
pronounced when re-running the test using only the nodes 
{0,2,4,6}. These were a set of four sensors at height 67 
cm closest to the bed, which achieve a RMS-median of 
0.83 when using all C = 5 channels. The improvement 
for this set, when increasing from one to four channels, is 
dramatic. The average RMS-median with C = 1 for this 
set of nodes is 2.54 - thus five channels achieves a 73% 
reduction. 

Interestingly, in some cases, a lower number of channels 
and a lower number of nodes can actually improve perfor- 
mance. We do not believe that it is better, on average, to 
use fewer sensors or fewer channels, because we may not 
know ahead of time which sensor positions will be best, 
and we can't know ahead of time which set of channels 
is best. Regardless, future work may exploit adaptive 
methods to select channels from among those that can be 
measured, or adapt which set of sensors should be operat- 
ing, in order to best estimate breathing rate. Such adaptive 
protocols could dramatically reduce the energy used in a 
breathing monitoring RF sensor network. Further, it would 
seem that better algorithms, which are more aligned with 
the statistics of the measured data, should be developed 
so that it doesn't degrade performance when more links' 
data is included. 



Fig. 10. RMS of the difference between tfie estimate and tfie sfiort- 
term median, as a function of S, the number of nodes, for nap 
experiment. 

To explore this further, we test specific (non-random) 
subsets of nodes and compare the results in Table|Vl First, 
we have eight nodes near the four corner of the bed, with 
nodes 0, 2, 4, and 6 at a height of about 0.67 m, and nodes 
1, 3, 5, 7 on the fioor. We first test the four nodes on the 
floor, {1, 3, 5, 7}, and see an RMS-median of 1.96 bpm. In 
contrast, the four nodes at 0.67 m, {0,2,4,6} have RMS- 
median of 0.72 bpm. We suspect that the links between 
nodes on the floor have signals which propagate largely 
along the floor, and as such do not interact strongly with 
the person in the bed. We also test using nodes {0, 1, 2, 3}, 
the nodes closest to the person. Using this set of nodes 
results in a very high RMS-median of 3.57 bpm. Although 
they are close to the person, they are close to the person's 
head, and no link line crosses the person's chest. 

B. Number of Channels 

In prior work |19) . nodes used only one channel. We 
are interested in whether our performance has improved 
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Number of Channels 

Fig. 11. RMS of the difference between the estimate and the short- 
term median, as a function of C, the number of channels, for the nap 
experiment. The results for using all twelve sensors {0, . . . , 11}, as 
well as using the node subset {0, 2, 4, 6}, are both shown. 
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VI. Conclusion 

This paper demonstrates key advances for the use of 
RSS in a wireless network to monitor the breathing of 
a person in the deployment area. First, we introduce a 
new method for mean removal that makes the method 
robust to small movements, a key requirement to be 
able to perform RSS-based breathing monitoring in real- 
world residential scenarios. Second, we show that the link 
PSD at the breathing rate can be used to estimate the 
location of the breathing person. We test both advances 
using a multi-channel network deployed in an apartment 
building. The results show the possibilities for accurate 
and reliable breathing monitoring and localization of a 
breathing person for a variety of applications. 

The area of wireless network breathing monitoring has 
many unanswered questions. What are good statistical and 
radio propagation models that explain a link's ability to 
measure breathing as a function of the person's position? 
Can adaptive algorithms and protocols be used to reduce 
the sampling requirements, and thus reduce the traffic 
and energy required for breathing monitoring? How should 
breathing be measured with other wireless hardware, such 
as 802.11 devices? What algorithms should be used to 
monitor and track the breathing of multiple people in the 
same deployment area? What WLAN protocols could be 
used to prevent an adversary from surreptitiously using 
someone's wireless network to eavesdrop on their breath- 
ing? Each of these questions may result in interesting and 
useful research directions. 
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